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We use the three dimensional uniformly frustrated XY model, as a model for a high temperature 
superconductor in an applied magnetic field, to explicitly measure the longitudinal correlation length 
£ z in the vortex line liquid phase. We determine the scaling of £ z with magnetic field and system 
anisotropy close to the vortex lattice melting transition. We apply our results to determine the 
extent of longitudinal correlations in YBCO just above melting. 
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It is now generally accepted that thermal fluctuations 
in the high T c superconductors lead, for a clean sample 
in the mixed state, to a first order melting of the vor- 
tex line lattice into a vortex line liquid. The properties 
of this vortex line liquid have been the subject of con- 
siderable investigation. An early theory by Feigel'man 
and co-workers jlj proposed that longitudinal (parallel 
to the applied field H ) superconducting coherence could 
still persist above melting. Flux transformer experiments 
in heavily twinned YBCO single crystals Q suggested 
support for this conclusion, as did early numerical simu- 
lations of the frustrated three dimensional (3D) XY 
model. However more recent experiments on untwinned 
YBCO single crystals by Lopez et al. ||] found longitu- 
dinal coherence to vanish simultaneously with melting, 
as have recent, more extensive, XY simulations by Hu et 
al. and by Nguyen and Sudb0 |6| . However it remains an 
important open question just how large the finite longitu- 
dinal correlations can become just above melting. Sim- 
ulations by Nordborg and Blatter within the "two 
dimensional (2D) boson" approximation, as well as gen- 
eral theoretical considerations predict a correlation 
length £ z ~ 7 _1 a v , where 7 = \ z /\±_ is the anisotropy 
ratio and a v = \J 4>o/ B is the average spacing between 
vortex lines. However analyses of experiments on un- 
twinned single crystal YBCO by Righi et al. and by 
Moore [[l0| have suggested that longitudinal correlations 
may be of the surprisingly larger micron scale. 

To investigate this issue, we carry out extensive new 
simulations of the frustrated 3D XY model for different 
values of applied flux density / and anisotropy n, ex- 
plicitly measuring the longitudinal correlation length £ z 
as determined by several different criteria. We find a 
good scaling of £ z with / and rj in the continuum limit, 
allowing us to estimate £, Z (T C ) in real YBCO single crys- 
tal samples. We find that longitudinal correlations at 
melting are enhanced with respect to the 2D boson ap- 
proximation, but not dramatically so. We also address 
several additional questions. We show, contrary to re- 
cent claims [ Q| , that there is only a single transition 
even in the isotropic model. In the very anisotropic limit 
£,z{T c ) < d, where a cross over to 2D behavior has been 



predicted |,]T§, we find no qualitative differences from 
the less anisotropic cases. We find that thermally excited 
vortex loops, which become important at low magnetic 
fields, can be described by an effective renormalization of 
the interaction between field induced vortex lines, and we 
find no evidence for a recently proposed transition within 
the vortex line liquid phase |i"3| , |l4| ]. 

Our model is the uniformly frustrated 3D XY model 
p5| , given by the Hamiltonian 



(1) 



1 . fi. 



where the sum is over the sites i of a cubic grid of points 
with unit basis vectors jl = x,y, z. 9i is the phase an- 
gle of the superconducting wavefunction on site i, and 
= (2tt/4>o) A • dl is the integral of the magnetic 
vector potential on the specified bond. The unit of the 
grid spacing along z is taken as d, the spacing between the 
weakly coupled CuO planes; the unit of the grid spacing 
in the xy plane is taken as £j_o, the bare vortex core size 
in the plane. The Hamiltonian (|l]) results from making 
the London approximation to the discretized Ginzburg- 
Landau energy, and assuming X/a v — > 00 so that the in- 
ternal magnetic field B can be taken as frozen and equal 
to the uniform applied field H. For a uniaxial anisotropic 
system with weak direction along z, the couplings are 
J x , y = J X = ^d/(16^ 3 Ai) and J z = ^l /(16n 3 X 2 z d), 
where Aj_ and X z are the penetration lengths in the re- 
spective directions. The anisotropy is given by the pa- 
rameter 
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and the magnetic field is taken uniform along z, with a 
density of flux quanta per plaquette of the grid, 
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/ and n are the two dimensionless parameters of our 
model. A more complete derivation of Eq. ([!]), and justi- 
fication for its use in modeling high T c materials, is given 
in Ref. B. Its advantage over the "2D boson" approxi- 
mation is in its more realistic vortex line interaction, and 
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in that it allows for the production of thermally activated 
vortex ring excitations, which may be important at small 
/■ 

To determine the relevant transitions in the model, we 
simulate Eq. ([!]) with periodic boundary conditions [jl6| , 
measuring the standard quantities [jl7|: (i) the helicity 
moduli parallel and perpendicular to the field, T z and 
Txj which measure phase coherence, and (ii) AS(K) = 
S(K) — S(R X [K]), where 5(kj.) is the average intraplanar 
vortex structure function, K is a reciprocal lattice vector 
of the ordered vortex lattice, and R x reflects K through 
the x axis; the difference is used so that AS vanishes in 
the liquid, and we average AS over the three smallest 
non-zero values of K. Our simulations at temperatures 
near the transition, for a lattice of size L± x L± x L z , 
consist typically of L\L Z Monte Carlo passes through 
the entire lattice for equilibration |l8|, followed by 2 x 
10 6 — 10 7 passes for computing averages. 

An example of our results is shown in Fig. |l| below for 
the case of isotropic couplings n = 1, and / = 1/20, for 
L± = 40 and several different sizes L z . If we denote the 
loss of longitudinal coherence, where Y z vanishes, as T c , 
and the melting of the vortex lattice, where AS vanishes, 
as T m , then only for the largest L z do we clearly observe a 
single transition with T m = T c jlji]. The strongest finite 
size effect is the increase in T m as L z increases. Our 
results explain recent simulations by Ryu and Stroud |ll]] 
which, using smaller systems, continued to suggest T m < 
T c for the isotropic model. Note that vanishes well 
below T c , indicating that the vortex lattice has depinned 
from our numerical grid well below its melting. 
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FIG. 1. Helicity moduli Y z , T z and vortex structure order 
parameter AS'(K) vs. T, for / = 1/20 and isotropic couplings 
r\ — 1, at system sizes L± — 40 and L z = 32, 64 and 128. Solid 
lines are guides to the eye. 

We next measure the longitudinal correlation lengths 
in the vortex line liquid, T c < T, as determined three 
different ways. The phase angle correlation length £ z 
and the vortex correlation length £ vz are defined by the 
correlation functions, 

C(z) = ^ e mr±,z)-e(r ± ,o)]) ^ e ~ z/ ^, T C <T (4) 
Cy(z) = (n z (r ± ,z)n z (r ± ,0)) ~ e~ 2 /«", T c < T (5) 



Here n z (r±, z) = ^[D x D6] • z is the vorticity in the xy 
plane at transverse position rj_ and height z (D is the lat- 
tice difference operator). We work in a gauge for which 
Ai z — 0. Our third length is determined by consider- 
ing the wavevector dependent helicity modulus Y z (kj_), 
which gives the linear response in supercurrent to a per- 
turbation in vector potential A z (k±) ]2C| ]. In D = 3, 
dimensional analysis gives T z ~ 1 /length. In the vortex 
liquid, provided one is not near any critical point where 
anomalous dimensions might come into play pl[ , T z (kjJ 
must vanish as k\ as k± — * 0. We therefore define the 
helicity correlation length £xz by 

r z (k ± ) = c£ Tz kl, T C <T (6) 

where c is a constant numerical factor, which we fix in 
an ad hoc manner by requiring £yz = £z at T = 1.2. 

For each case we have considered, we first carefully 
choose L z sufficiently large so as to observe a single sharp 
first order melting transition; however L z must not be too 
large, in order that we are still able to cool into the vor- 
tex lattice state without getting trapped in a supercooled 
liquid. For such a value of L z , we carefully monitor the 
time sequence of AS and determine T c as the tempera- 
ture at which the system seems to be switching equally 
between vortex lattice and vortex liquid states. To ac- 
curately measure correlation lengths, we then repeat the 
simulations with a larger value of L z £ Z (T C ), cooling 
down to the predetermined T c . 

In Fig. H we show results for £ 2 , £ vz and £x z vs. T, for 
the case of / = 1/20, n 2 = 9, with Lj_ = 40 and L z = 128 
(for T > 0.8, L z = 64). We also show the specific heat C, 
as computed from energy fluctuations. The peak in C at 
Tc2 n is identified as the cross-over where, upon cooling, 
local superconducting order first develops Q| . 




FIG. 2. Phase, vortex, and helicity correlation lengths, £ 2 , 
£ vz , and £y z vs. T for / = 1/20 and rj 2 — 9, for system size 
L± = 40 and L z = 128. Also shown is the specific heat C. 
Solid lines are guides to the eye. 

We see that all three lengths increase similarly as one 
cools towards T c . No noticeable feature is seen near T C 2. 
£ 2 is slightly larger than £ vz by a factor of about 1.3. 
The numerical factor of Eq. (^|) is found to be c = 74. In 
determining £ z and £ vz from Eqs. (^) and (||), we fit our 
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data self consistently within the range £ z < z < L z /3, 
averaging over the position rj_. In determining £yz, we 
fit Eq. (^) to the two smallest non-zero values of k± , since 
we found that T z (kx) quickly saturated to a constant 
as kj_ increased. This unfortunately limits the accuracy 
with which we can determine £y2- Some examples of 
our fits for £ z and £xz are shown in Fig. || below. In 
the remainder of this work we now focus on the phase 
correlation length £ z . 
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FIG. 3. a) Phase correlation C(z) vs. z, and b) helicity 
T z (k^) vs. k±, for several different values of T for the pa- 
rameters of Fig. p. The solid lines are fits to Eqs. (Q) and (g) 
that determine £ z and (t z . 



We therefore expect that the dimensionless £ z /d should 
scale as 1/d, and so we conclude £ z /d ~ 1/^/Jrf. In 
Fig. H below we replot the results of Fig. f|as d) y/ frf 
vs. T/T c . We see that as T — ► T c , most of the data 
collapse to a single curve. Deviations from this curve 
represent situations when either £ z /d is small ~ O(l), or 
when T is sufficiently large (approaching T C 2) that ther- 
mally excited vortex rings start to dominate the total 
vorticity of the system. In the first case, the discrete- 
ness of our grid spacing d clearly becomes an important 
length scale. In the second the density of ther- 

mal rings is determined by the vortex core energy, and 
hence by the core sizes d and £j_o, again the discreteness 
of our grid becomes evident. Such deviations thus occur 
for all cases at sufficiently high T, and also for the case 
/ = 1/32, r; 2 = 1250 at all T. This latter case was specifi- 
cally chosen so that, according to continuum expressions, 
one would expect £ Z (T C ) < d and so to be in the so-called 
"2D" limit of very weakly coupled layers |p|,^2|. We see 
that for this case £ z (T c )/d lies below the other data, indi- 
cating an even smaller correlation length than one would 
expect in a continuum. However we otherwise found no 
anomalous behavior for this case: there remained only a 
single transition where longitudinal coherence and vortex 
lattice order vanished simultaneously. 



In Fig. |i| we show our results for £ z vs. T, for several 
different values of the parameters / and rj. 
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FIG. 4. Phase correlation length £ 2 vs. T for parame- 
ter values / = 1/20 and rj 2 — 1, 9, 40 (system sizes are 
L± = 40 and L z = 192, 128, 32 respectively); / = 1/12 and 
rj 2 = 9 (system size is L± — 24 and L z — 64); / = 1/32 
and rj 2 = 1250 (system size is L± = 64 and L z — 8); 
and / = 1/100 and rj 2 — 1 (system size is L± — 100 and 
L z = 256). 

We can now argue as follows for the dependence of £ z 
on / and rj. For small values of /, such that £xo «C a v , 
we expect that £ z should be independent of the vortex 
core size £j_o. From Eqs. (§-|J) we see that the only com- 
bination of / and r\ that is independent of £j_o is frf 1 . 
Furthermore, for large £ z we expect our discretizing grid 
to become a reasonable approximation of the continuum 
and so £ z should be independent of the layer spacing d. 
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FIG. 5. Data of Fig. |replotted as 



T/T c . 



Except for the "2D" case discussed above, our other 
data, when appropriately scaled as in Fig. ^ all coincide 
at T c . We therefore conclude that, for these cases, our 
model is well approximating continuum behavior. From 
the specific numerical value of (£ z /d)y/ fi] 2 at T c in Fig. ^ 
we can therefore conclude that in a uniaxial anisotropic 
superconductor in the "3D" continuum limit, 



5.57 a v 



(7) 



where the second equality follows from Eqs. (0-|3j), with 
7 = \ z /\± the anisotropy, and a v the average spacing 
between vortex lines. Applying this result to YBCO, 
for which 7 ~ 7, we conclude that £ Z (T C ) ~ 0.86a v , or 
£ z ~ 0.023/i for a field of B ~ H = 4T. 

Our result above may be compared with that of re- 
cent "2D boson" simulations of Nordborg and Blatter 
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M, which yielded £, VZ (T C ) — 1.7"/~ 1 a v . If we take from 
Fig. H that £z{T c ) — l.3£ vz (T c ), then the more realistic 
vortex line interaction of the XY model gives roughly a 
2.5 fold increase in £ 2 (T C ) over the boson model. This 
remains, however, well below the micron scale. 

Except for the case / = 1/100, our numerical re- 
sults are all in the limit of sufficiently large /, such 
that T c (f, 77) lies well below the zero field critical point 
T c (0,?7). For these cases, therefore, thermally excited 
vortex rings are not playing any significant role at our 
melting transitions. One way to see this is to note that 
for these cases, the melting temperatures T c (/, rj) obey 
quite well the expectation of the Lindemann criterion 
(which ignores thermal rings), T c /J± cx Jrf B. To 
see this, note our result above that ^ z (T c )/d oc l/\Zfrf, 
and hence, if the Lindemann criterion holds, we expect 
[£ z (T c )/d\/[T c /J±] to be a constant. In Fig. g we see 
that the loci of points ([£ 2 (T c )/d], [T c /J±]) do indeed lie 
on quite close to a straight line intersecting the origin. 

The case / = 1/100, however, clearly lies off this line. 
This is as expected: as / decreases and T c (f, 77) increases, 
one eventually enters the critical region of the / = tran- 
sition, where thermally excited rings play a significant 
role in renormalizing the effective interactions between 
the magnetic field induced vortex lines, and suppress 
the melting transition below the value predicted by the 
"bare" Lindemann criterion, so that limy^o T c (f, rj) — 
T c (0, 77). In this case, our argument that £z(T c ) should 
be independent of £j_o becomes less obvious. Neverthe- 
less, we see in Fig. || that the data for / = 1/100 agrees 
quite well with our scaling assumption in the near vicin- 
ity of T c . We therefore conclude that the main effect 
of the thermally excited rings at melting is indeed ad- 
equately described by a renormalization of vortex line 
couplings [p2| , and so our result of Eq. (Q) will continue 
to hold in the low field region, although with a possible 
renormalization of the anisotropy parameter 7. 

Recently, Tesanovic Q has argued that there may 
still be a singular vortex ring blowup transition at low 
fields, within the normal vortex line liquid. Recent XY 
simulations by Nguyen and Sudb0 [^4) have claimed to 
identify this transition in terms of a sharp percolation 
transition of transverse vortex paths, which takes place 
in the vicinity of the peak in the specific heat. They 
find that, within the XY model, this percolation transi- 
tion is more clearly distinct from the melting transition 
at larger rather than smaller fields. One of Tesanovic's 
predictions is that the length £yz will have a discontinu- 
ous decrease at this transition. This prediction has been 
one of our motivations in computing £yz- In Fig. || we 
find no clear evidence for such behavior in £xz- It there- 
fore remains unclear, if such a percolation or ring blowup 
transition does exists, whether it has any noticeable effect 
on thermodynamically measurable quantities. 
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